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ABSTRACT 


i'  This  report  presents  the  findings  and  conclusions  of  a  computer  simu¬ 
lation  study  of  the  Welland  Canal  and  proposed  alternatives  to  the  Welland. 
The  objective  of  this  study  was  to  establish  the  limits  of  service  for 
alternative  structural  and  nonstructural  improvements  in  terms  of  delay 
and  transit  time.  Waterborne  transport  demand  through  the  year  2030  was 
considered  in  the  analysis,  demand  being  represented  by  two  factors, 
traffic  density  and  fleet  composition.  The  results  of  the  study  were  por¬ 
trayed  as  a  series  of  transit  time  response  curves  for  each  configuration, 
plotted  as  a  function  of  the  transport  demand. 
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FOREWORD 

The  work  described  in  this  interim  report  was  performed  by  the  Penn¬ 
sylvania  Transportation  and  Traffic  Safety  Center  (PTTSC)  at  The 
Pennsylvania  State  University  for  the  U.S.  Army  Corps  of  Engineers,  North 
Central  Division,  under  contract  number  DACW-23-72-C-0066.  The  contract 
period  is  from  1  July  1972  to  31  June  1973. 

The  effort  reported  herein  represents  a  continuation  of  work  at  PTTSC 
in  the  general  area  of  waterway  systems  planning  and  analysis.  An  earlier 
interim  report  entitled  Great  Lakes  Simulation  Studies,  Volume  1 — NETSIM: 

A  General  Network  Simulator  documented  the  development  of  a  simulation 
model  for  multiple  channel  deep  draft  navigation  studies.  The  present 
report  applies  this  model  to  the  Lake  Erie-Lake  Ontario  navigation  subsys¬ 
tem. 

PTTSC  personnel  contributing  to  this  investigation  were  the  following: 

Joseph  L.  Carroll,  Principal  Investigator 

Head,  Transportation  Systems  Division,  PTTSC  and 

Professor  of  Business  Administration 

Srikanth  Rao 

Research  Assistant  in  Management 

Michael  S.  Bronzini 

Research  Assistant  in  Civil  Engineering 

David  C.  Nowading 

Research  Assistant  in  Business  Logistics 

John  E.  McLauchlan 

Graduate  Assistant  in  Business  Logistics 


Richard  Dunhoff 

Graduate  Assistant  in  Business  Logistics 
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I.  INTRODUCTION 

This  report  is  part  of  a  three-phased  study  commissioned  by  the  U.S. 
Corps  of  Engineers  to  assist  them  in  their  appraisal  of  the  capacity  of  the 
Great  Lakes  navigation  system  with  focus  on  the  structural  and  nonstruc- 
tural  improvements  that  may  be  postulated  by  the  increasing  waterway  demand. 
The  first  two  tasks  were  concerned  with  the  capacity  of  the  existing 
Welland  Canal  and  the  potential  need  for  a  parallel  Lake  Erie-Lake  Ontario 

(LE-LO)  Canal  via  the  Niagara  River.  The  final  task  was  formulated  to 

assist  the  Corps  of  Engineers  with  additional  Great  Lakes  System  Studies 
including  St.  Lawrence  River  Navigation,  Great  Lakes  Season  Extension,  and 
Sault  Ste  Marie  (Sabin)  locks.  These  tasks  may  be  summarized  as  follows: 

1.  to  develop  a  LE-LO  Navigation  Simulation  model, 

2.  to  apply  the  model  to  simulation  studies  of  the  Welland  Canal  and 

proposed  alternatives  to  the  Welland, 

3.  to  revise  the  initial  simulation  model  so  as  to  Include  the 
capabilities  needed  for  comprehensive  Great  Lakes-St.  Lawrence 
System  simulations. 

Task  1,  the  development  of  the  simulation  model,  has  been  documented  pre¬ 
viously  (1).*  The  present  report  presents  the  findings  of  Task  2,  viz., 

2 

the  application  of  this  model  to  the  Welland-Niagara  complex. 

A.  TASK  2:  LE-LO  -  WELLAND  CANAL  STUDIES 

This  phase  of  the  study  involves  four  major  subtasks,  as  follows: 


^Numbers  in  parentheses  refer  to  the  list  of  References. 
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The  terms  "Niagara  Canal"  and  "LE-LO  Canal"  are  used  interchangeably  in 
this  report. 
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1.  to  establish  the  expected  limits  of  service  of  the  existing 
Welland  Canal, 

2.  to  establish  the  expected  incremental  increase  in  service  poten¬ 
tial  of  the  existing  Welland  Canal  under  assumptions  of  improved 
locking  procedures  and  an  improved  traffic  control  system, 

3.  to  determine  the  expected  performance  of  a  combined  Welland- 
Niagara  system  with  configuration  alternatives  of  four,  five  and 
six  locks  in  series  in  the  Niagara  Canal  in  combination  with  the 
existing  Welland  Canal, 

4.  to  examine  the  expected  performance  of  a  replacement  for  the 
Welland  Canal  consisting  of  a  series  of  four  super  locks  plus  a 
guard  lock  towards  the  mouth  of  Lake  Erie. 

These  configurations  were  subjected  to  current  and  anticipated  levels  of 
traffic,  fleet  composition,  ship  size,  and  operating  procedures.  The 
primary  measure  of  system  performance  was  system  transit  time.  This  vari¬ 
able  reflects  both  the  service  levels  provided  by  system  facilities  and 
any  delays  that  occur  due  to  congestion.  In  addition,  measures  of  lock 
utilization,  lock  processing  time,  and  time  spent  in  queues  were  taken  so 
that  the  system  response  could  be  stated  in  terms  of  delays  due  to  conges¬ 
tion  and  lock  utilization.  However,  no  analysis  of  the  effects  of  delays 
and  system  congestion  upon  demand  was  undertaken.  Hence,  the  emphasis  of 
the  study  was  placed  upon  determining  what  configurations  of  navigation 
facilities  are  required  to  meet  the  prospective  transportation  demand  and 
enable  the  network  to  function  effectively  as  a  system. 
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B.  SYSTEM  DESIGN  ALTERNATIVES 

Simulation  of  the  existing  Welland  Canal  established  the  calibration 
values  for  the  model's  parameters.  Recent  traffic  data  for  the  Welland 
Canal  were  used  to  establish  state-dependent  relationships  between  vessel 
transit  time  and  the  number  of  ships  in  the  canal.  The  canal  was  simu¬ 
lated  using  input  data  representative  of  the  conditions  corresponding  to 
those  for  which  the  traffic  data  were  compiled.  The  capacity  of  the  exist¬ 
ing  canal  was  established  by  subjecting  it  to  waterway  transport  demand 
through  the  year  2030. 

It  has  been  reported  that  improvements  in  locking  operations  and  the 

installation  of  a  traffic  control  system  in  the  Welland  Canal  have  led  to 

an  increase  greater  than  33  percent  in  the  potential  number  of  lockages 

per  day  (2).  The  study  cited  also  states  that  on  the  average,  1  minute 

3 

saved  per  lock  cycle  saves  1  hour  in  round-trip  transit  time  for  the  ves¬ 
sels.  An  earlier  report  stated  that  a  lock  cycle  time  of  70  minutes  might 
eventually  be  achieved  (3).  The  conclusions  reached  in  these  two  reports 
suggest  that  the  already  improved  locking  time  might  be  further  reduced  by 
an  additional  2  minutes  resulting  in  a  further  improvement  in  round-trip 
transit  time  on  the  order  of  2  hours. 

The  effects  of  nonstructural  improvements  in  the  Welland  Canal  can 
be  inferred  under  the  assumption  that  the  optimization  of  vessel  scheduling 
and  locking  procedures  can,  in  fact,  lead  to  improvements  in  vessel  transit 
time  on  the  order  of  those  stated  above.  Simulation  runs  for  future  time 
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A  lock  cycle  is  the  time  needed  for  a  lock  to  move  one  vessel  down  and 
one  vessel  up.  This  includes  all  the  ship  movements  and  lock  operation 
from  the  time  a  boat  is  instructed  to  enter  a  lock  to  the  time  when  the 
lock  is  ready  to  receive  the  next  lockage  in  the  same  direction. 
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periods  using  inputs  revised  to  incorporate  these  efficiencies  were  made 
in  order  to  ascertain  the  capacity  of  an  "improved"  Welland  Canal. 

For  the  purposes  of  this  study,  the  Welland  Canal  was  modeled  as  a 
set  of  six  entities,  corresponding  with  the  traffic  data  that  were  compiled, 
where  the  operations  within  each  entity  were  inferred  rather  than  specifi¬ 
cally  modeled.  Transit  time  through  the  system  depended  on  the  number  of 
ships  in  the  canal  and  was  calibrated  on  the  basis  of  empirical  data.  The 
effects  of  additional  nonstructural  improvements  to  the  canal  were  deter¬ 
mined  by  simulation  experiments  as  described  above. 

Four  different  structural  alternatives  for  the  Welland  Canal  were 
studied.  Three  of  these  alternatives  consisted  of  two  parallel  channels 
in  a  combined  Welland-Niagara  system.  The  performance  of  this  system  was 
measured  for  configurations  of  four,  five,  and  six  locks  in  series  in  the 
Niagara  Canal  in  combination  with  the  existing  eight-lock  Welland.  The 
simulation  used  transit  time  data  representative  of  the  future  improved 
operations  at  the  Welland  and  performance  data  for  the  Niagara  based  on 
specifications  for  locks  and  reaches  provided  by  the  Corps. 

The  fourth  structural  configuration  involves  a  unilateral  Canadian 

alternative  to  the  Welland  where  the  existing  system  would  be  replaced  by 

4 

a  series  of  five  locks  of  greater  lift  and  1200'  x  100'  dimensions.  This 
configuration  was  simulated  as  a  single  channel  using  performance  data  for 
locks  and  reaches  from  specifications  provided  by  the  Corps. 

C.  METHODOLOGY 

The  basic  methodology  for  the  Welland-Canal  simulation  experiments 
entailed  the  division  of  traffic  between  parallel  facilities.  This  factor 


4 

Includes  a  guard  lock. 
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dictated  that  a  channel  assignment  mechanism  be  incorporated  in  the 
modeling  system.  As  a  result  of  discussions  with  the  Corps,  the  groundwork 
for  this  system  was  laid  in  a  previous  study  conducted  at  PTTSC  (4).  The 
following  model  specification  was  formulated: 

1.  Canal  Operating  Rule:  all  ships  may  use  either  existing  or  new 
facilities  when  physically  possible;  ships  too  large  to  use  the 
existing  locks  must  be  assigned  to  the  large  new  locks. 

2.  Lock  Operating  Rules:  ships  queuing  on  both  sides  of  a  lock 
should  be  serviced  alternatively.  A  recycle  "lookahead"  capabil¬ 
ity  should  be  included  such  that,  there  being  no  queues  at  a  lock, 
the  lock  water  level  should  be  adjusted  to  accept  the  second  of 
two  ships  traveling  in  the  same  direction,  prior  to  its  arrival 

at  the  lock,  time  permitting. 

3.  Reach  Operating  Rule:  ships  should  be  allowed  to  catch  up  with, 
or  fall  further  behind  a  preceding  ship  but  should  not  be  allowed 
to  pass  a  preceding  ship  in  a  canal  reach. 

4.  Assignment  Decision  Rule:  ships  should  be  assigned  between  paral¬ 
lel  facilities  on  the  basis  of  the  least  expected  transit  time. 

The  model  formulation  has  been  extensively  documented  in  the  Task  1  report 
but  a  brief  mention  of  the  assignment  decision  mechanism  is  made  here  since 
it  has  a  significant  bearing  on  the  methodology  of  the  simulation  experi¬ 
ments. 

The  assignment  decision  involves  the  simulation  of  each  parallel 
branch  separately,  using  a  special  option  embedded  in  the  simulation  model 
to  produce  an  "experience  data  base"  (EDB).  These  experience  data  bases 
are  then  statistically  analyzed  to  select  the  most  significant  variables 
and  to  establish  coefficient  values  for  use  in  a  set  of  transit  time 
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estimating  functions  for  each  branch.  The  coefficient  values  and  variable 
identities  for  use  in  each  branch  estimating  equation  are  specified  by  the 
user.  With  these  calibrated  functions  available,  simulation  of  the  sys¬ 
tem  of  parallel  canals  can  be  performed.  When  a  ship  arrives  at  a  channel 
choice  point,  the  channel  assignment  is  based  on  the  least  expected  transit 
time. 

This  empirical  approach  was  adopted  because  it  was  found  that  an 
analytical  approach  to  the  determination  of  expected  transit  time  through 
a  multiple  lock  and  reach  canal  was  extremely  complex  and  was  intractable. 
It  is,  therefore,  postulated  that  a  statistical  relationship  should  exist 
between  the  conditions  existing  in  a  canal  when  a  ship  arrives  at  the 
assignment  decision  point  and  the  time  that  will  be  required  to  subsequent¬ 
ly  travel  through  that  canal.  By  performing  a  simulation  (called  an  EDB 
run)  of  a  given  channel  configuration  (i.e,,  a  specific  canal  branch),  it 
is  possible  to  develop  an  experience  data  base  which  includes  prior  canal 
conditions  for  each  ship  arrival  and  subsequent  ship  transit  time.  This 
data  base  may  then  be  analyzed,  using  a  standard  statistical  program,  to 
establish  the  required  relationships  between  expected  transit  time  and 
canal  conditions.  Separate  relationships  for  each  of  the  Welland-Niagara 
configurations,  differentiated  by  direction  of  travel,  have  been  developed 
and  are  presented  in  Appendix  B. 

Thus,  the  EDB  concept  requiring  preliminary  simulation  runs  for  the 
formulation  of  the  expected  transit  time  functions  constitutes  the  first 
phase  of  the  methodology.  The  actual  simulation  of  the  network  configura¬ 
tions,  using  these  expected  transit  time  functions,  for  various  input 
factors  through  different  levels  of  transport  demand  forms  the  second 
phase.  An  actual  simulation  (called  an  EVENT  LOG  run)  generates  a  log  of 
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each  event  occurrence  during  simulation,  and  this  event  log  is  usually 
placed  on  an  external  output  device  such  as  a  magnetic  tape  for  permanent 
storage.  It  is  this  event  log  which  is  used  during  the  third  phase  of  the 
methodology . 

Because  of  the  complexity  and  sheer  volume  of  the  simulation  model's 
logic  structures,  the  burden  of  statistical  evaluation  was  shifted  to  a 
post-simulation  phase.  Thus,  the  third  phase  consists  of  an  event  log 
post-processor  which  utilizes  the  event  log, generated  during  the  second 
phase  to  summarize  significant  statistical  data.  The  task  of  this  event 
log  processor  is  conceptually  simple  since  it  merely  reads  in  the  event 
log  as  input  data  and  translates  it  into  periodic  (thereby  providing 
intermediate  snapshots)  statistical  summaries.  The  processor,  however, 
also  includes  options  for  determining  such  simulation  specifics  as  tran¬ 
sient  time  and  autocorrelation  in  data. 

The  distinction  between  "transient  state"  and  "steady  state"  is  ex¬ 
tremely  relevant  for  simulation  studies.  A  transient  state  is  defined  as 
a  condition  in  which  operating  characteristics  of  a  system  are  changing 
with  time  in  irregular  fashion;  and,  conversely,  a  steady  state  means  an 
environment  in  which  system  operating  characteristics  are  not  time  vary¬ 
ing.  In  the  simulation  context,  a  transient  state  usually  exists  at  the 
start  of  a  simulation  run  when  the  system  is  in  some  initial  state.  During 
the  course  of  the  simulation,  the  system  eventually  reaches  a  steady  state 
as  it  accepts  more  entities,  although  under  high  utilization  the  system 
may  depart  from  steady  state. 

Thus,  the  amount  of  simulation  time  required  to  reach  this  steady 
state  (called  warm-up  time)  is  of  critical  importance  to  the  analyst  as  he 
will  normally  want  to  make  decisions  based  only  on  steady  state  performance 
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measures.  The  warm-up  time  may  be  determined  in  a  number  of  ways  (5). 

The  event  log  processor  uses  a  procedure  that  examines  the  behavior  of 
delays  at  locks  through  spectral  analysis. 

The  spectral  analysis  methodology  was  originally  necessitated  by  the 
presence  of  autocorrelation  in  simulation  data.  That  is,  in  most  simula¬ 
tions  generating  time  series  data  there  is  a  certain  degree  of  dependence 
on  past  events  which  violates  a  common  statistical  requirement  that  the 
data  be  independent.  A  spectral  analysis  approach,  however,  can  be  used 
in  such  cases  to  measure  autocorrelation  and  take  it  into  account  in  the 
subsequent  analysis.  This  approach  has  been  documented  elsewhere  (6); 
hence,  no  description  is  given  here.  The  use  of  this  approach  in  the 
Welland-Niagara  studies  led  to  the  choice  of  5000  minutes  as  the  warm-up 
time  and  30,000  minutes  as  the  total  simulation  time. 

D.  ASSUMPTIONS 

All  simulation  experiments  were  conducted  under  certain  assumptions 
agreed  to  between  PTTSC  and  the  Corps  of  Engineers.  These  assumptions  are: 

1.  The  system  never  breaks  down. 

2.  All  locks  operate  under  the  "S0QA"  (Serve  Opposing  Queues  Alterna¬ 
tively)  rule,  where  ships  queuing  on  both  sides  of  a  lock  will  be 
served  alternately,  reverting  to  "First  Come  First  Served"  only 
when  one  queue  becomes  empty. 

3.  The  channel-choice  decision  rule  for  parallel  branches  is  based 
on  least  expected  transit  time. 

4.  There  are  no  double,  tandem  or  combination  lockages. 

5.  There  is  no  priority  service  given  to  any  ship. 

6.  Time  of  day  does  not  affect  traffic  levels  (nor  day  of  week,  nor 
month  of  season). 
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7.  Locks  operate  24  hours  per  day,  7  days  per  week. 

8.  There  will  be  no  passing  in  any  reach. 

9.  Reach  transit  time  will  not  be  a  function  of  ship  size  or  ship 
direction  of  travel. 

10.  Statistics  are  gathered  under  steady  state  conditions. 

11.  Arrivals  are  random  at  the  endpoints  of  the  canal  in  Lake  Erie 
and  Lake  Ontario. 
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II.  INPUT  DATA  BASE 


A.  DATA  SOURCES 

The  data  source  for  the  Welland  Canal  was  the  St.  Lawrence  Seaway 
Authority  which  supplied  "Welland  Canal  Vessel  Transit  Analysis  Daily 
Details"  for  the  months  of  April,  June,  August,  and  October,  1971.  The 
data  for  the  months  of  April  and  October  were  influenced  by  seasonality 
effects  but  preliminary  examinations  indicated  that  August  data  would 
adequately  serve  the  purposes  of  this  study. 

The  data  sources  for  the  Niagara  Canal  were  the  Corps  of  Engineers, 
Buffalo  and  North  Central  districts.  The  data  provided  by  the  Corps  con¬ 
sisted  of  system  configuration  parameters,  transit  time  distributions  for 
locks  and  reaches,  and  finally  the  fleet  data  for  the  Welland-Niagara 
studies. 

B.  DATA  DESCRIPTION 

A  complete  description  of  the  data  used  in  this  study  is  shown  in 
Appendix  D.  This  section  merely  provides  some  of  the  salient  features  of 
the  input  data. 

1.  Fleet  Data 

The  fleet  data  consisted  of  actual  average  traffic  level  in  1971 
through  the  Welland  Canal  and  projected  average  levels  for  the  years  1980, 
1985,  1990,  1995,  2000,  2010,  2020,  and  2030.  Two  estimates  for  each  of 
the  projected  levels  were  given,  one  representing  normal  growth  and  the 
other  representing  accelerated  growth  as  shown  in  Table  1. 

Associated  with  this  traffic  projection  was  a  fleet  composition  fac¬ 
tor  representing  the  estimated  trend  towards  larger  fleet  size.  Fleet 
data  consisted  of  three  vessel  classes,  delineated  by  length  of  vessel  as 


follows: 


11 


TABLE 


1 .  AVERAGE 

DAILY  VESSEL 

TRANSITS  (ships/day) 

Normal 

Growth 

Accelerated 

Growth 

1970 

25.50 

25.50 

1980 

26.50 

27.40 

1985 

27.00 

28.55 

1990 

27.50 

29.70 

1995 

28.00 

31.15 

2000 

28.50 

32.60 

2010 

29.50 

35.70 

2020 

30.90 

39.80 

2030 

32.20 

44.30 

Note:  Data  for  1970  are  actual,  for  the  others  are 
projected. 


* 
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Class  I  :  1  -  399  feet 

Class  II  :  400  -  730  feet 

Class  III:  731  -  1150  feet 

A  description  of  the  projected  fleet  composition  (percentage  distribution 

by  class)  for  the  years  1970  through  2030  Is  given  In  Table  2.  It  may 
be  noted  that  Class  III  Is  comprised  of  vessels  too  large  to  be  processed 
through  the  Welland  Canal. 

2 .  Lock  Data 

The  simulation  model,  hereon  referred  to  as  NETSIM/SHIP  (NETwork 
SIMulator  of  SHIP  movements),^  represents  a  lock  operation  in  terms  of 
nine  elements  as  depicted  In  Figure  1  and  delineated  in  Table  3.  Each  of 
these  nine  elements  can  be  described  in  the  simulation  model  by  either  a 
frequency  distribution  or  an  average  time  of  operation.  The  data  for 
these  elements  by  class  size  and  direction  of  travel  is  shown  in  Table  4. 
All  the  locks  in  the  Niagara  Canal  of  four-,  five-,  and  six-lock 
configurations  as  well  as  the  four-super-lock  Welland  configuration  were 
of  size  1200'  x  110'  and  used  these  identical  data. 

3.  Welland  Transit  Data 

The  only  major  treatment  of  raw  data  occurred  for  the  "Welland  Canal 
Vessel  Transit  Analysis  Daily  Details"  for  the  month  of  August.  This 
treatment  had  two  objectives. 

The  first  objective  was  to  obtain  a  transit  time  distribution  for 
each  of  the  Canal's  six  entitles  for  calibration  purposes.  Figure  2 
provides  a  schematic  description  of  the  six  entities  representing  the 


^A  complete  description  of  the  simulation  model  is  given  in  Reference  (1). 
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TABLE  2.  FLEET  COMPOSITION-PERCENTAGE  DISTRIBUTION  BY  CLASS 


Class  I 
(l'-399') 

Class  II 
(400'-730') 

Class  III 
(731'-1150’ 

1970 

30.00 

70.00 

0.00 

1980 

23.30 

74.70 

2.00 

1985 

19.85 

74.60 

5.55 

1990 

16.40 

74.50 

9.10 

1995 

14.00 

73.45 

12.55 

2000 

11.60 

72.40 

16.00 

2010 

9.60 

68.60 

21.80 

2020 

6.60 

65.00 

28.40 

2030 

5.30 

60.00 

34.70 

Note:  Data  for  1970  are  actual,  for  the  others  are  projected 
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Figure  1.  Schematic  of  Simulated  Locking  Time  Events 
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TABLE  3.  DESCRIPTION  OF  LOCKING  TIME  EVENTS 


EVENT  NAME 


DESCRIPTION 


A.  MOVING  ENTRY 

Begins: 

Ends: 

B.  QUEUED  ENTRY 

Begins : 
Ends. 


into  the  lock  chamber  from  the  Clear  Point  at 
the  und  of  the  entry  throat. 

when  the  bow  of  the  ship  passes  the  Clear  Point, 
when  the  gates  begin  to  close  astern  of  the 
ship. 

into  the  lock  chamber  from  the  head  of  the 
queue  adjacent  to  the  Clear  Point, 
when  the  gates  are  fully  open  and  the  chamber, 
is  free. 

when  the  gates  begin  to  close  astern  of  the 
ship. 


C.  MOVING  APPROACH 
ENTRY 

Begins : 
Ends: 


from  the  Clear  Point  to  a  position  in  the  entry 
throat  just  clear  of  the  entry  gate, 
when  the  bow  of  the  ship  passes  the  Clear  Point, 
when  the  ship  comes  to  rest  in  the  entry  throat. 


D.  STATIONARY  APPROACH 
ENTRY 


Begins : 
Ends: 


from  the  head  of  the  queue  adjacent  to  the  Clear 
Point  to  a  position  in  the  entry  throat  just 
clear  of  the  entry  gate. 

when  the  bow  of  the  ship  passes  the  Clear  Point, 
when  the  ship  comes  to  rest  in  the  entry  throat. 


E.  SHORT  ENTRY 


Begins : 
Ends: 


into  the  lock  chamber  from  a  stationary  position 
just  clear  of  the  entry  gate  in  the  entry 
throat. 

when  the  gates  are  fully  open  and  the  chamber 
is  free. 

when  the  gates  begin  to  close  astern  of  the  ship. 


F.  LOCKAGE 

Begins : 

Ends: 

G.  CHAMBER  EXIT 

Begins: 
Ends : 


of  a  ship  at  rest  in  the  chamber. 

when  the  entry  gates  begin  to  close  astern  of 

the  ship. 

when  the  exit  gates  are  fully  opened  after  the 
change  in  water  level. 

from  chamber  to  position  where  the  ship's  stern 
clears  the  exit  gate. 

when  the  exit  gates  are  fully  opened  after  the 
change  in  water  level. 

when  the  ship's  stern  is  clear  of  the  exit  gate. 
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TABLE  3.  CONTINUED 


DESCRIPTION 

from  position  where  the  ship's  stern  clears  the 
exit  gate  to  the  Clear  Point  at  the  end  of  the 
exit  throat. 

when  the  ship’s  stern  is  clear  of  the  exit  gate, 
when  the  stern  of  the  ship  passes  the  Clear 
Point. 

I.  RECYCLE  of  the  water  level  with  no  ship  in  the  chamber. 

Begins:  when  the  gates  begin  to  close. 

Ends:  when  the  opposite  gates  are  fully  open  to  re¬ 

ceive  an  incoming  ship. 


EVENT  NAME 
H.  THROAT  EXIT 

Begins: 

Ends: 
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TABLE  4.  LOCK  DATA  FOR  WELLAND-NIAGARA  STUDIES 
(Average  Times  of  Operation  in  Minutes) 


Lock  Elements 

Direction 

Class  I 

Class  II 

Class  III 

A. 

MOVING  ENTRY 

Up 

21 

28 

37 

Down 

20 

26 

34 

B. 

QUEUED  ENTRY 

Up 

24 

31 

40 

Down 

21 

29 

38 

C. 

MOVING  APPROACH 

16 

20 

ENTRY 

Up 

14 

Down 

12 

14 

19 

D. 

STATIONARY 

19 

23 

APPROACH  ENTRY 

Up 

17 

Down 

15 

17 

22 

E. 

SHORT  ENTRY 

Up 

12 

16 

23 

Down 

11 

15 

23 

F. 

LOCKAGE  (PROCESS) 

Up 

14 

16 

18 

Down 

14 

16 

18 

G. 

CHAMBER  EXIT 

Up 

8 

9 

10 

Down 

7 

8 

10 

H. 

THROAT  EXIT 

Up 

6 

7 

8 

Down 

5 

7 

8 

I. 

RECYCLE 

Up 

11 

11 

11 

Down 

11 

11 

11 

19 


Welland  Canal,  and  Table  5  summarizes  the  corresponding  transit  data. 

The  second  objective  was  the  derivation  of  a  transit  time  equation  as 
a  function  of  system  conditions  which  would  be  used  to  predict  the  ex¬ 
pected  transit  time  through  the  canal  during  simulation.  The  results  of 
this  treatment  are  given  in  Appendix  B.  Note  that  the  corresponding 
transit  time  equation  for  the  proposed  Niagara  Canal  had  to  be  derived 
through  an  EDB  simulation  run;  the  results  of  which  are  also  given  in 
Appendix  B.  An  EDB  run  was  not  necessary  for  the  Welland  Canal,  however, 
since  empirical  data  were  available. 


from  Lock  1  to  Lock  7 
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III.  CALIBRATION 


A.  THE  WELLAND  CANAL 

The  simulation  run  for  the  existing  Welland  Canal  network,  modeled  as 
a  set  of  six  reaches  under  a  1971  traffic  load,  served  as  the  base  run  for 
subsequent  analysis  of  various  other  configurations.  This  simulation 
methodology  did  not  use  Monte  Carlo  sampling  from  random  probability 
distributions  for  a  vessel's  transit  time,  but  rather  utilized  the  empiri¬ 
cal  transit  time-system  condition  relationship  obtained  through  regression 
as  explained  in  Appendix  B. 

Thus,  at  the  moment  of  a  vessel's  entrance  into  the  Welland  Canal, 
its  expected  canal  transit  time  was  computed  as  a  function  of  the  existing 
system  conditions.  Transit  through  individual  reaches  was  subsequently 
accomplished  through  random  sampling  from  prior  distributions  whose  means 
were  empirically  determined  fractions  of  the  expected  canal  transit  time. 

Since  the  technique  described  above  constitutes,  in  principle,  a 
projection  of  the  empirical  relationship  rather  than  true  Monte  Carlo 
simulation,  it  will  be  hereon  referred  to  as  "ETT  simulation."  The 
projection  of  this  relationship  is  valid  for  only  finite  deviations  from 
the  1971  traffic  load,  in  fact  only  for  the  range  of  the  empirical  data 
from  which  the  relationship  had  been  obtained.  Since  the  traffic  load 
for  some  of  the  future  periods  in  the  analysis  falls  outside  this  range, 
queuing  theory  as  explained  in  Appendix  C  was  used  to  supplement  ETT 
simulation. 

The  calibration  run,  then,  consists  of  an  ETT  simulation  of  the 
existing  Welland  Canal  under  a  1971  traffic  load.  A  fleet  mix  of  30  percent 
Class  I  and  70  percent  Class  II  was  used.  Results  of  the  simulation  run, 
shown  in  Table  6,  are  differentiated  by  direction  and  are  given  in  both 
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minutes  and  hours.  Transit  times  through  individual  sections  of  the  canal 
compare  favorably  with  the  empirical  data  (see  Table  5).  The  latter  are 
understated  by  about  4  percent  for  downstream  transit  and  by  about  1.5  percent 
for  upstream  travel  with  an  overall  underestimation  of  about  2  percent. 

B.  THE  PROPOSED  NIAGARA  CANAL 

In  order  to  determine  the  relationship  between  the  Welland  Canal 
and  the  proposed  Niagara  Canal,  the  minimum  time  to  transit  the  latter  for 
its  various  configurations  was  calculated  from  the  input  data  as  follows: 


Transit 
through  Locks 
(minutes) 

Transit 

through  Reaches 
(minutes) 

Total 

(minutes) 

UP  DOWN 

UP 

DOWN 

4-Lock  Niagara 

240  228 

527 

767 

755 

5-Lock  Niagara 

300  285 

506 

806 

791 

6-Lock  Niagara 

360  342 

481 

841 

823 

These  figures  assume  an  average  locking  time  of  60  minutes  including 
entry  and  exit  and,  in  fact,  represent  the  input  data  for  Class  II  vessels. 
The  minimum  times  shown  in  this  table  assume  that  all  locks  would  be 
properly  set  to  receive  a  vessel  without  any  delay.  Under  such  circum¬ 
stances,  the  minimum  transit  time  through  the  Niagara  Canal  for  all 
configurations  (except  upstream  transit  for  the  four  lock  Niagara)  is  still 
greater  than  the  mean  transit  time  through  the  Welland  Canal  under  normal 
traffic.  This  suggests  that  in  the  subsequent  simulations,  the  proposed 
Niagara  Canal  might  serve  as  the  backup  alternative  for  a  Welland  under 
congestion,  except  for  Class  III  vessels  which  would  be  too  large  to  have 
a  channel  choice.  This  is  precisely  borne  out  in  the  later  studies 
although  traffic  tends  to  equalize  in  the  long  run  due  to  the  changes  in 


fleet  composition. 
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IV.  THE  SINGLE  WELLAND  STUDIES 

Under  chis  classification,  there  are  again  three  configurations: 

1.  the  existing  Welland  with  no  further  major  improvements, 

2.  the  improved  Welland  with  nonstructural  Improvements, 

3.  the  improved  Welland  with  structural  improvements  consisting 
of  four  new  super  locks  of  1200'  x  110'  dimensions  and, 

in  addition,  a  guard  lock  to  compensate  for  fluctuations 
with  Lake  Erie.^ 

Since  the  study  of  each  network  was  unique  in  itself,  no  general  intro¬ 
duction  is  supplied  except  to  state  the  primary  need  of  these  studies  was 
to  establish  the  expected  limits  of  service  of  each  network. 

A.  ESTABLISHING  THE  CAPACITY  OF  THE  EXISTING  WELLAND 

A  reasonable  estimate  of  the  capacity  of  the  Welland  Canal  may  be 
derived  from  a  number  of  ways.  Three  techniques  that  prominently  come 
to  mind  are  Monte  Carlo  simulation,  analytical  methods  such  as  queuing 
theory  and  regression  analysis,  and  a  qualitative  appraisal  based  on  the 
nominal  capacities  of  individual  locks  in  the  Welland. 

The  utility  of  Monte  Carlo  simulation  models  for  waterway  systems 
analysis  has  been  previously  demonstrated  (? ) .  The  special  case  of  a 
sophisticated  control  system  at  Welland  poses  problems  for  such  simula¬ 
tion,  however.  In  particular,  rules  and  procedures  such  as  the  priority 
rating  for  ships  proceeding  toward  lock  7  over  those  which  have  passed  it, 
the  preferential  treatment  to  faster  ships  if  a  gap  develops  in  the 
traffic  toward  lock  7,  special  rules  for  heavy  one-way  traffic  and  for 
ships  most  susceptible  to  wind  and  fog  are  difficult  to  model. 


^Existing  locks  would  be  leveled  (Source:  St.  Lawrence  Seaway  Authority). 
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For  these  and  other  reasons,  Monte  Carlo  simulation  has  not  been  employed 
in  the  Welland  studies. 

Discussions  with  the  St.  Lawrence  River  Authority  have  indicated  that 
a  maximum  of  40  lockages  per  day  can  be  accomplished  at  the  Welland  locks. 
This  accomplishment  can  only  occur  under  conditions  of  favorable  weather, 
of  large  queues  (in  the  order  of  10)  on  both  sides,  and  of  normal 
operation.  Any  deviations  from  these  corv it ’ ons  such  as  a  fog  or  an 
accident  leads  to  a  rapid  increase  in  the  queue.  A  more  reasonable  level 
of  capacity  quoted  is  around  36  lockages  per  day  or  about  90  percent  of 
nominal  capacity.  Increasing  the  capacity  m  terms  of  vessels  per  day 
or  lockages  per  day  is,  however,  only  one  measure  of  performance.  What 
matters  most  to  shipping  companies  is  the  total  canal  travel  time,  including 
both  the  time  they  wait  to  enter  the  canal  and  the  actual  canal  transit 
time  itself.  This  total  travel  time  could  be  worsened  even  though 
capacity  in  terms  of  lockages  per  day  is  increased.  For  example,  lock 
utilizations  of  the  order  given  above  require  large  queues  on  either  side 
of  the  lock,  and  hence  it  is  necessary  to  have  many  vessels  inside  the 
canal  at  any  one  time.  Then,  even  Lhough  the  input  rate  and  output  rate 
from  the  canal  were  increased,  the  canal  transit  times  associated  with 
this  level  of  operation  could  also  be  extremely  large. 

The  capacity  of  the  Welland  Canal  may  also  be  estimated  by  way  of 
analytical  methods.  Two  such  methods  mentioned  above  and  employed 
in  this  study  are  the  queuing  and  regression  (ETT  simulation)  models. 

Results 

2 

Given  the  arrival  rate  in  terms  of  expected  lockages  per  day 

2 

This  characterizes  one  of  the  assumptions  of  this  study,  viz.,  there 
will  be  no  tandem  or  multiple  lockages.  In  reality,  these  constituted  about 
11  percent  of  the  arrival  rate  in  1971,  a  decline  from  20  percent  in  1965. 


26 

and  the  largest  lock  service  time,  the  ratio  of  these  two  parameters 
provides  a  measure  of  system  capacity.  When  this  ratio  reaches  one 
or  100  percent,  the  system  is  said  to  have  found  its  nominal  capacity. 

Discussions  with  the  St.  Lawrence  River  Authority  revealed  that 
lock  7,  which  is  the  bottleneck  on  the  Welland  Canal,  had  a  lock  cycle 
of  72  minutes  or  a  single  lockage  time  (including  ship  entry  and  exit) 
of  36  minutes.  Solving  for  the  arrival  rate  in  terms  of  lockages  per 
day  needed  to  reach  the  nominal  system  capacity,  this  figure  turns  out 
to  be  exactly  40  (lockages  per  day),  no  doubt  the  source  of  the  earlier 
information. 

In  a  similar  manner  the  expected  system  capacity  in  percentage  for 
each  year  of  projected  growth  can  be  calculated,  and  these  results  are 
shown  in  Table  7  ror  both  normal  and  accelerated  rates.  These  figures 
are  lower  than  the  lock  utilizatioris  that  are  actually  realized  because 
the  effects  of  coordinating  tandem  and  multiple  lockages  are  not  included. 

If  the  latter  continue  to  develop  in  future  traffic,  the  canal  could 
saturate  earlier  than  Table  7  would  indicate. 

3 

The  results  of  the  ETT  simulation  runs  are  shewn  in  Table  A1  and 
additionally  are  plotted  along  with  the  results  ot  the  queuing  model  in 
Figure  3.  In  observing  the  ETT  curve,  the  statement  concerning  the  range 
of  the  empirical  data  over  which  the  regression  equation  was  fitted  is 
once  more  brought  to  light. 

To  draw  inferences  regarding  Welland  capacity,  it  is  necessary  to 
set  forth  a  definition  of  capacity.  (It  is  assumed  here  that  the  definition 
that  is  of  interest  would  consider  the  transit  time  associated  with  some 
particular  level  of  input  traffic  or  lock  utilization  rather  than  the 
nominal  capacity  that  has  been  referred  to  previously.)  One  can  obtain 

3 

Collected  in  Appendix  A. 


TABLE  7.  THE  EXISTING  WELLAND  CANAL  -  CAPACITY  SUMMARY 
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YEAR 

Figure  3.  Transit  Time  Measures  for  the  Existing  Welland  Canal 
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some  assistance  on  this  point  from  the  use  of  75  percent  of  nominal  capacity 
in  literature  as  an  indication  of  an  economical  or  practical  level  ol 
capacity  (for  example,  see  (8),  (9),  (10)).  If  this  criterion  is  used  to 
indicate  capacity,  it  is  met  under  accelerated  growth  projections  as  early 
as  1990  and  under  normal  growth  projections  at  2010. 

In  Figure  3,  these  two  measures  fall  in  those  portions  of  the  curves 
where  the  slopes  are  changing  rapidly.  This  indicates  that  further  traffic 
can  be  introduced  into  the  canal;  however,  the  delays  associated  with 
such  additional  increases  would  rapidly  accelerate.  Further,  the  traffic 
figures  used  for  1990  accelerated  and  2010  normal  projections  are  29.7  and 
29.5  lockages  per  day,  respectively,  and  represent  the  seasonal  daily 
average  for  the  Welland  Canal.  Traffic  during  certain  periods  of  the  season 
may  be  greater  than  these  seasonal  averages,  however. 

The  statement  "Welland  capacity  is  reached  in  1990  under  accelerated 
projections  and  in  2010  under  normal  growth  projections,"  therefore, 
involves  the  following  assumptions: 

1.  Capacity  used  in  the  above  statement  is  75  percent  of  the 
maximum  theoretical  capacity. 

2.  The  traffic  projections  used  are  seasonal  averages. 

3.  The  traffic  projections  used  are  the  number  of  lockages  per  day 
and  not  vessels  per  day. 

4.  The  effects  of  coordinating  tandem  and  multiple  lockages  are 
not  Included. 

5.  Structural  as  well  as  nonstructural  operations  as  in 
existence  in  1971  are  assumed  to  remain  unchanged. 

The  assumptions  given  above  also  apply  to  the  next  two  studies. 
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B.  ESTABLISHING  THE  CAPACITY  OF  THE  IMPROVED  WELLAND  (NONSTRUCTURAL) 

In  the  ETT  simulation  of  this  task,  the  effects  of  various  non- 
structural  improvements  were  assumed  to  have  resulted  in  first,  a  reduction 
of  the  lock  cycle  from  72  minutes  to  70  minutes;  and  second,  the  transla¬ 
tion  of  this  efficiency  into  a  net  transportation  savings  of  2  hours  in 
round-trip  transit  time. 

Thus,  the  simulation  itself  consisted  of  incorporating  these 
efficiencies  in  the  ETT  model's  transit  time  distributions  derived  for 
each  vessel  when  it  arrives  at  the  call-in-point  and  therefore,  is  again 
an  extrapolation  of  a  current  empirical  relationship  to  future  demand. 
Hence,  this  technique  suffers  from  the  same  shortcoming  as  the  previous 
one  and  is  valid  only  for  finite  deviations  from  current  capacity. 

Capacity  is  related  to  the  input  stream  as  shown  in  Table  8.  The 
capacity  figures  represent  a  ratio  of  the  input  rate  to  the  lock  service 
rate  where  the  latter  assumes  a  lock  service  time  of  35  minutes  per 
lockages.  Nominal  capacity  for  this  system  is  41  lockages  per  day,  but 
using  the  criterion  from  the  previous  section,  75  percent  of  system 
capacity  is  reached  at  a  vessel  arrival  rate  of  about  31  lockages  per  day. 
An  estimate  of  the  transit  time  associated  with  these  levels  of  system 
utilization  can  be  obtained  from  Figure  4.  This  figure  compares  the 
simulation  and  queuing  model  results;  the  similarity  of  behavior  between 
this  and  the  preceding  network  is  as  expected. 

C.  ESTABLISHING  THE  CAPACITY  OF  THE  IMPROVED  WELLAND  (STRUCTURAL) 

A  preliminary  study  of  this  network  (for  a  diagram,  see  Appendix  D) 
consisted  of  a  Monte  Carlo  simulation  using  the  locking  data  provided  by 
the  Corps.  The  network  was  subjected  to  traffic  demands  for  years  1980, 
1985,  and  1990.  Further  experimentation  was  deemed  unnecessary  as  the 
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Figure  4.  Transit  Time  Measures  for  the  Nonstructurally 
Improved  Welland  Canal 
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system  was  judged  to  have  reached  its  capacity.  The  primary  reason  for 
this  was  the  data  providing  for  an  average  lock  cycle  of  120  minutes  (or 
a  single  lockage  time  of  60  minutes  including  vessel  entry  and  exit) . 

At  this  service  rate,  nominal  system  capacity  is  only  about  24  lockages 
per  day. 

Subsequent  discussions  with  the  St.  Lawrence  Seaway  Authority 
indicated  that  a  lock  cycle  time  of  80  minutes  might  well  be  attainable. 

The  queuing  model  was  therefore  employed  to  develop  a  series  of  response 
curves  where  different  lock  cycle  times  were  employed.  The  results  of 
this  analysis  are  shown  in  Figure  5,  where  the  different  curves  represent 
lock  cycle  times  ranging  from  72  minutes  to  96  minutes. 

These  results  stress  the  ramarkable  sensitivity  of  system  capacity 
to  lock  cycle  times  and  suggests  that  significant  economics  exist  in 
reductions  in  say  lock  approach  and  entry  time.  Large  vessel  sizes 
permit  an  increase  in  tonnage  capacity;  however,  the  resultant  loss  of 
maneuverability  affecting  a  more  cautious  and  slower  lock  approach  leads 
to  a  decrease  in  system  capacity.  This  tradeoff  indicates  that  signi¬ 
ficant  implications  for  fleet  and  facility  design  may  be  obtained  from 
the  type  of  analysis  made  above.  The  simulation  model  is  particularly 
useful  in  this  context  since  it  contains  a  detailed  representation  of  the 
lock  and  can,  therefore,  provide  fast  computations  of  the  various  tradeoffs 
as  mentioned  above — assuming  of  course  that  the  relationship  between 
vessel  size  and  locking  time  can  be  quantified  properly  as  input  data. 


» 


AVERAGE  TRANSIT  TIME  (Minutes) 


Figure  5.  Transit  Performance  for  the  Structurally 
Improved  Welland  Canal 
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V.  TWIN  CANAL  STUDIES 


A.  SIMULATION  RESULTS 

Three  configurations  were  simulated  in  these  studies.  All  configura¬ 
tions  had  the  Welland  Canal  in  common,  but  the  Niagara  Canal  consisted  of 
either  four,  five, or  six  locks  in  series.  The  latter,  therefore,  also 
contained  either  seven,  eight,  or  nine  reaches  as  shown  in  the  network 
diagrams  in  Appendix  D. 

Passing  was  not  permitted  in  any  of  the  reaches  except  the  end 
reaches  of  the  Niagara  Canal  which  form  the  outliers  of  the  lakes.  Vessels 
originate  at  ports  and  arrive  at  channel  assignment  decision  nodes  where 
NETSIM/SHIP's  channel  assignment  mechanism  is  set  in  motion  to  determine 
the  channel  offering  the  least  expected  transit  time. 

The  simulation  experiments  subjected  each  network  to  increasing 
transport  demand  from  1980  through  to  2030,  if  necessary,  in  five  year 
increments  up  to  year  2000  and  in  ten  year  increments  thereafter.  Each 
simulation  was  examined  for  signs  of  saturation  to  determine  if  the 
next  higher  level  of  demand  experiment  was  necessary. 

The  principal  parameters  of  interest  in  these  simulations  were 
transit  time,  delays,  and  lock  utilizations.  The  latter  two  measures 
could  only  be  observed  at  the  locks  in  the  Niagara  Canal,  and  performance 
through  the  Welland  Canal  was  assessed  in  terms  of  transit  time. 

1.  Welland-Four  Lock  Niagara  Subsystem 

The  overall  behavior  of  this  system  can  be  captured  in  a  graphical 
display  of  average  transit  time  as  shown  in  Figure  6.  Certain  other 
selected  system  parameters  are  exhibited  in  Table  A. 2. 


The  tendency  toward  curvature  in  Figure  6  may  be  inferred  from  a 
closer  examination  of  system  parameters  for  each  canal.  The  number  of 
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Figure  6.  Average  Transit  Time  for  the  Welland- 
Four  Lock  Niagara  Subsystem 
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journeys  and  the  average  transit  time  through  each  canal  are  given  in 
Table  A. 3.  The  average  transit  time  through  the  Welland  drops  from  the 
calibration  run  value  of  757  minutes  (12  hours  and  36  minutes)  in  1971  to 
673  minutes  (11  hours  and  13  minutes)  in  1980.  Throughout  the  rest  of  the 
period  of  experimentation,  the  average  transit  time  levels  off  because  the 
decreasing  traffic  through  Welland  falls  outside  the  range  of  the  empiri¬ 
cal  data  from  which  the  ETT  relationship  had  been  obtained.  One  can  note 
however,  that  in  conjunction  with  the  decrease  in  the  Welland  traffic,  the 
composition  of  the  fleet  becomes  dominated  in  time  by  Class  II  vessels. 

The  load  factor  through  the  Niagara  Canal,  on  the  other  hand,  increases 
significantly  from  period  to  period,  this  being  in  large  part  due  to 
the  projected  increase  in  vessels  of  Class  III. 

The  system  performance  measures  for  the  Niagara  Canal  in  terms  of 
delays  and  lock  utilization  are  given  in  Table  9.  The  canal  is  extremely 
underutilized  and  is  probably  not  in  steady  state  until  about  1990.  In 
year  2030  with  about  12  ships  per  day,  average  canal  delay  is  about 
2-1/2  hours  with  lock  utilizations  approximately  60  percent. 

2.  Welland-Five  and  Six  Lock  Niagara  Subsystems 

The  simulation  results  for  these  configurations  were  very  similar 
to  those  for  the  Welland-Four  Lock  Niagara  network.  The  main  items  of 
departure  from  the  latter  were  the  degree  of  traffic  diversion  from 
the  Welland  to  the  Niagara  in  Classes  I  and  II  and  th''  average  transit 
time  through  the  Niagara  Canal. 

The  average  system  transit  times  for  these  configurations  are 
plotted  in  Figure  7.  The  occurrence  of  high  delays  at  high  utilization 
rates,  contributing  to  the  curvature,  is  again  established  in  Tables  10 
and  11  and  in  Tables  A. 4,  A. 5,  A. 6,  and  A. 7. 
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TABLE  9.  THE  NIAGARA  CANAL  PERFORMANCE  MEASURES  FOR  THE 

WELLAND-FOUR  LOCK  NIAGARA  SUBSYSTEM  UNDER  NORMAL  GROWTH 


Year 

Average  Canal 
Delay 
(min) 

Total  Canal 
Delay 
(min) 

Average  Lock 
Utilization 
(%) 

Maximum  Delay 
at  any 
Single  Lock 
(min) 

1980 

2.8 

76 

6.4 

47 

1985 

9.7 

397 

10.5 

65 

1990 

9.2 

452 

12.8 

95 

1995 

18.1 

1,226 

18.1 

101 

2000 

36.9 

3,539 

24.4 

134 

2010 

47.2 

5,941 

29.2 

251 

2020 

84.6 

13,802 

45.6 

305 

2030 

141.7 

29,598 

55.4 

427 

r 
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Figure  7.  Average  Transit  Time  for  the  Welland-Five  and 
Six  Lock  Niagara  Subsystems 
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The  behavior  of  the  five-lock  Niagara  usually  lay  In  between  those 
of  the  four-  and  six-lock  configurations  although  no  statistical  difference 
could  be  established  in  many  cases. 

B.  DISCUSSION  OF  RESULTS 

The  results  of  the  Monte  Carlo  simulations  suggest  that  any  one  of 
the  twin-canal  configurations  would  function  adequately  in  meeting 
transport  demand  up  to  and  including  that  for  2030.  Although  the 
Welland-Four  Lock  Niagara  network  seems  to  offer  lower  transit  time  than 
the  others,  no  absolute  statistical  difference  between  the  performances 
could  be  established.  Therefore,  the  desirability  of  a  particular 
configuration  must  rest  on  other  economic  criteria  such  as  initial 
capital  outlay  and  operating  and  maintenance  expenditures. 

If  the  system  interactions  remain  the  same,  and  if  the  locking  time 
remains  constant,  one  can  speculate  on  the  saturation  arrival  rates, 
i.e.,  that  level  of  traffic  which  would  give  rise  to  high  lock  utiliza¬ 
tions.  Given  the  assumptions  above,  the  lock  utilizations  can  be 
calculated  as  the  ratio  of  the  arrival  rate  to  the  lock  service  rate  and 
as  such  are  given  below. 


Lock  Utilization 

Mean  Service  Time  1/y 

Arrival  Rate  A 

75% 

72  minutes 

15.0  ships /day 

80% 

72 

15.99 

90% 

72 

18.0 

99% 

72 

19.8 

The  above  relationships  are  exact  and  indicate  that  if  19.8  ships 
per  day  arrive  into  the  Niagara  Canal  (not  the  Welland-Niagara  subsystem) 
where  the  projected  average  lock  service  time  is  72  minutes,  then  the 
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locks  will  experience  a  99  percent  utilization  rate.  To  obtain  the  life 

K 

of  a  Welland-Niagara  twin-canal  subsystem,  it  is  further  necessary  to 
consider:  1)  the  capacity  of  the  Welland,  2)  potential  increases  in 
capacity  through  further  nonstructural  improvements,  and  3)  trends  in 
fleet  composition. 

In  summary,  the  simulations  suggest  that  the  primary  benefits  from 
the  construction  of  a  twin  canal  are  transportation  savings  not  due  to 
large  reductions  in  total  canal  transit  time,  but  rather  due  to  reduced 
ship  idle  time  in  combination  with  the  increased  capacity  permitted  by 
larger  lock  dimensions. 


43 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  SUMMARY 

The  following  points  summarize  the  conclusions  from  the  individual 
network  simulations.  A  pictorial  summary  is  given  in  Figures  8  and  9. 

1.  Performance  statistics  from  simulations  of  the  twin-canal  networks 
suggested  that  any  of  the  four-,  five-,  or  six-lock  Niagara  Canal  in 
combination  with  Welland  would  perform  equally  well  under  projected 
traffic  up  to  year  2030.  This  conclusion  must  be  tempered  by  the 
following  observations: 

i.  All  twin-canal  simulations  used  a  single  set  of  locking  data. 

Under  these  conditions,  no  statistical  differences  could  be 
established  between  any  of  the  three  configurations.  Changes 
in  locking  data  among  the  various  networks  could  revise  this 
result. 

ii.  In  connection  with  the  above  and  specifically  with  regard  to 
the  alternative  Niagara  Canal  configurations,  the  traffic  was 
completely  dominated  by  Class  III  vessels  and  as  such  an 
average  lock  processing  time  of  70  minutes  was  obtained. 

Bearing  in  mind  that  this  figure  refers  to  a  1200'  x  110' 
lock  and  that  the  lock  data  were  based  on  original  observa¬ 
tions  at  Poe  Lock,  comparison  may  be  made  with  the  following 
lock  operations: 

Identification  Dimension  Average  Lift  Average  Time  Per  Lockage 

existing  locks  860'  x  80'  43 '-48'  36  minutes 

Welland  with  nonstructural  860'  x  80'  43'-48'  35  min.  (expected) 

Canal  improvements 

super  locks  1200'  x  110'  80'  40  min.  (expected) 

Eisenhower-Snell  860'  x  80'  38'-49'  37  min.  (for  large 

vessels) 


Figure  8.  Transit  Summary  under  Normal  Growth 


Figure  9.  Transit  Summary  under  Accelerated  Growth 
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Thus,  any  major  improvement  in  lock  design  (or  operations) 
invalidating  the  use  of  the  70-minute  locking  time  could  give 
rise  to  significant  increases  in  benefits  in  the  form  of  reduced 
delay  and  higher  capacity. 

iii.  The  simulation  runs  assumed  a  uniform  and  smooth  operation  in 
the  channels  (reaches)  connecting  the  locks.  No  channel 
delays  were  incorporated  into  the  model.  This  is  not  a  grave 
problem  for  the  twin-canal  studies  since  the  systems  were  never 
in  serious  congestion.  At  higher  utilization  rates,  however, 
major  channel  delays  could  "deregulate"  vessel  interarrival 
patterns  leading  to  larger  delays, 
iv.  The  study  assumed  a  SOQA^  (Serve  Opposing  Queues  Alternately) 
rule  and  further  included  a  service  lookahead  feature  so  that 
no  "unnecessary"  ship  delays  occurred. 

2.  Given  the  locking  data  for  the  Niagara  Canal,  immediate  benefits  of 
a  twin-channel  relative  to  the  calibration  run  (viz.,  the  existing 
Welland)  were  manifested  in  the  form  of  reduced  delays  and  increased 
capacity  but  not  in  reduced  canal  transit  time.  This  is  readily 
evident  when  one  considers  the  minimum  time  required  to  transit  the 
Niagara  under  the  null  condition  as  earlier  shown  on  page  23. 

3.  Given  the  input  data  for  the  Niagara  Canal,  the  ETT  mechanism 
channeled  the  bulk  of  the  traffic  in  the  early  years  of  traffic 
projection  into  the  Welland  Canal  with  the  Niagara  Canal  serving 
primarily  as  a  back  up  channel.  Traffic  through  the  latter  consisted 
primarily  of  Class  III  vessels,  and  delays  constituted  a  minor 
portion  of  the  total  transit  time.  Indeed,  the  locks  reached  50  percent 


Previous  studies  at  PTTSC  have  demonstrated  that  under  a  variety  of  system 
conditions  the  SOQA  rule  leads  to  a  significantly  improved  waterway  opera¬ 
tion  than  a  FIFO  (First  In  First  Out)  rule.  For  example  see  (11). 
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utilizations  in  2020  under  accelerated  growth  and  only  in  year  2030 
under  normal  growth. 

4.  Queuing  models,  when  calibrated,  proved  quite  useful  in  performing 
sensitivity  analysis,  i.e.,  in  testing  fluctuations  in  system 
parameters.  Selected  systems,  particularly  in  the  Single  Welland 
studies,  were  shown  to  be  extremely  sensitive  to  the  magnitude  of 
lock  cycles. 

5.  The  term  nominal  capacity  refers  to  the  maximum  theoretical  capacity 
and  is  determined  by  the  maximum  number  of  vessels  that  can  be 
locked  through  in  any  time  period.  The  occurrence  of  random  arrivals, 
however,  imposes  an  increasing  delay  cost  as  the  system  approaches 
nominal  capacity.  Thus,  what  is  needed  is  a  practical  capacity  which 
takes  into  account  this  delay  function.  The  definition  of  practical 
capacity  is  to  a  certain  extent  arbitrary  and  is  truly  a  function  of 

a  number  of  factors,  each  of  which  may  be  different  for  various 
individuals  in  various  systems  under  various  circumstances.  For  the 
purposes  of  this  study,  practical  capacity  was  defined  as  that  point 
when  the  system  reaches  75  percent  of  its  nominal  capacity. 

6.  Nominal  capacity  of  the  existing  Welland  was  shown  to  be  40  lockages 
per  day.  Practical  capacity  as  defined  above  was  reached  as  early  as 
1990  under  accelerated  growth  and  as  late  as  2010  under  normal  growth. 

7.  Nonstructural  improvements  in  the  Welland  Canal  leading  to  a  reduction 
in  lock  cycle  significantly  increased  the  economic  capacity  of  the 
system.  Given  the  traffic  projections,  the  capacity  of  the  system  was 
increased  by  almost  a  decade  under  normal  growth  from  that  given  above. 

8.  For  the  Welland  Super  Locks  subsystem,  a  lock  cycle  of  about  120  minutes 
as  given  in  the  input  data  was  determined  to  be  unreasonably  high. 
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Consequently,  a  set  of  canal  transit  curves  for  varying  lock  service 
rates  was  developed.  The  capacity  of  this  system  was  found  to  be 
extremely  sensitive  to  lock  service  times. 

9.  All  the  studies  assumed  the  absence  of  tandem  and  multiple  lockages. 
Although  these  currently  constitute  only  10-15  percent  of  the  total 
in  the  Welland  Canal,  the  relaxation  of  this  assumption  could  lead 
to  higher  lock  utilizations  than  those  predicted  in  the  studies. 

10.  An  overall  survey  of  the  analysis  suggests  that  nonstructural 

improvements,  particularly  those  related  to  lock  efficiencies  (and 
smooth  channel  operation) ,  need  to  be  given  more  emphasis  in  light 
of  their  impact  on  system  capacity. 

B.  EVALUATION 

Any  evaluation  of  the  individual  studies  and  th'  ‘.r  results  as  described 
above  must  a  priori  recognize  the  limitations  and  assumptions  of  the  studies. 
The  intent  of  the  following  presentation  is  to  make  the  reader  aware  of 
these  limiting  factors  so  that  he  may  assess  the  quality  and  correctness 
of  the  studies. 

A  primary  limitation  of  the  study  is  that  it  ignores  the  interaction 
between  system  congestion  and  transport  demand.  To  overcome  this  limita¬ 
tion,  it  would  be  necessary  to  chart  the  sensitivity  of  individual 
commodity  traffic  to  system  performance,  in  effect  forming  a  feedback 
loop  structure  that  allows  a  study  of  the  economics  of  alternative  facili¬ 
ties.  Lacking  this  overall  evaluation  model,  two  levels  of  transport 
demand  were  used  for  each  year  of  future  forecast,  thus  providing  what 
is  expected  to  be  the  range  of  actual  demand. 

The  simulations  performed  in  this  study  are  nondynamic  in  the  sense 
that  the  nature  of  system  Interaction  as  programmed  into  the  existing 
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models  Is  not  time-dependent.  For  example,  the  simulation  model  does  not 
forecast  the  future  status  of  the  system.  Future  conditions  may  be 
simulated  only  by  forecasting  the  future  values  of  the  simulation  inputs, 
as  vp s  done  for  this  study. 

A  third  limitation  concerns  the  precision  of  results.  Given  that 

the  main  objective  of  the  study  was  to  establish  system  performance  for 

alternative  levels  of  transport  demand,  and  given  the  constraints  of 

time  and  resources  available,  no  further  attempt  to  estimate  precise 

2 

delay  values  was  made.  Intermediate  simulation  results  or  other  means 

3 

of  replication  are  required  to  obtain  statistically  sound  delay  estimates. 
Treating  these  simulation  results  as  precise  numerical  values  is  tantamount 
to  ignoring  the  stochastic  features  of  both  the  models  and  the  system. 

Finally,  the  simulation  results  are  conditional  upon  the  accuracy 
of  traffic  forecast.  Transport  demand  for  future  years  were  prepared  by 
NCD  under  the  assumption  that  both  average  vessel  size  and  frequency  of 
trips  would  gradually  increase.  While  this  procedure  is  certainly  reason¬ 
able,  a  different  assumption  or  different  growth  rates  would  undoubtedly 
produce  different  results. 

Within  the  limitations  discussed  above  and  in  light  of  the  conclusions 
from  the  individual  subsystem  studies,  the  following  evaluation  may  be 
expounded . 


2 

In  connection  with  this,  one  may  note  that  delay  values  were  obtainable 
only  for  the  Niagara  Canal  studies. 


3 

Note  that  the  simulation  methodology  provides  for  a  permanent  record, 
so  that  simulations  need  not  be  duplicated. 
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The  most  clear-cut  result  of  these  studies  is  that  twin-canal  con¬ 
figurations  are  able  to  provide  better  service  over  a  longer  period  of  time 
than  the  single  Welland  configurations.  The  Welland-Niagara  subsystems 
distinctly  reflect  excess  capacity  through  the  end  of  the  current 
millennium.  The  benefits  in  terms  of  reduced  delays  and  the  absence 
of  congestion  may  have  been  understated,  however,  if  there  exists  a 
potential  for  significant  improvements  in  lock  cycle  times. 

The  single  Welland  studies  provide  testimony  to  the  volatile  relation¬ 
ship  between  system  capacity  and  the  lock  cycles.  While  the  structurally 
improved  Welland  (the  Welland  Super  Locks  subsystem)  affords  some  relief 
in  the  long  run  under  the  assumptions  of  the  analysis,  it  too  merits 
further  attention  by  way  of  1)  effect  of  enlarged  lock  dimensions  on 
trip  frequencies  and  2)  the  potential  improvement  in  transit-time 
through  traffic  control.  Neither  of  these  parameters  could  be  explicitly 
incorporated  into  the  models  and  must  be  either  qualitatively  assessed  or 
expressed  in  terms  of  efficiencies  in  model  input  data. 

The  simulation  results  presented  in  this  report  intimate  the  need 
for  some  additional  analysis  within  the  framework  of  an  overall  evaluation 
of  the  Welland  Canal.  This  analysis  should,  in  particular,  verify  the 
accuracy  of  the  simulation  data  and  explore  the  sensitivity  of  system 
performance  to  incremental  changes  in  facility  service  rates.  These 
changes  in  turn  may  elucidate  the  implications  for  fleet  and  facility 
designs  and  for  system  operation  procedures.  Given  this  type  of  sensi¬ 
tivity  analysis,  the  potential  benefits  of  improvements  may  be  more  readily 
assessed. 
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APPENDIX  A 

SELECTED  SIMULATION  RESULTS 


2030  5.3  94.7  16,021  32.3  1,060  44.3  1,217 


54 


TABLE  A. 2.  SELECTED  SYSTEM  PARAMETERS  FOR  THE  WELLAND-FOUR 
LOCK  NIAGARA  SUBSYSTEM  UNDER  NORMAL  GROWTH 


Vessel  Journeys 

Total  Transit  Time 
(minutes) 

Average  Transit  Time 
(minutes) 

1980 

458 

313,715 

684.9 

1985 

470 

328,953 

699.9 

1990 

475 

333,925 

703.4 

1995 

480 

342,528 

712.7 

2000 

491 

356,318 

725.0 

2010 

508 

381,508 

751.5 

2020 

527 

411,060 

780.6 

2030 

554 

454,834 
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TABLE  A. 4.  SELECTED  SYSTEM  PARAMETERS  FOR  THE  WELLAND-SIX 
LOCK  NIAGARA  SUBSYSTEM  UNDER  NORMAL  GROWTH 


Vessel 

Journeys 

Total  Transit  Time 
(minutes) 

Average  Transit  Time 
(minutes) 

1980 

458 

322,234 

703.5 

1985 

470 

334,810 

712.3 

1990 

475 

341,373 

718.6 

1995 

480 

353,968 

737.4 

2000 

491 

367,120 

747.7 

2010 

508 

393,308 

774.2 

2020 

527 

428,411 

812.9 

2030 

554 

474,073 

855.7 
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TABLE  A. 6.  SELECTED  SYSTEM  PARAMETERS  FOR  THE  WELLAND-FIVE 
LOCK  NIAGARA  SUBSYSTEM  UNDER  NORMAL  GROWTH 


Vessel 

Journeys 

Total  Transit  Time 
(minutes) 

Average  Transit  Time 
(minutes) 

1980 

458 

314,646 

687 

1985 

470 

329,047 

700 

1990 

475 

334,875 

705 

1995 

480 

343,680 

716 

2000 

489 

360,393 

737 

2010 

510 

385,050 

755 

2020 

528 

419,760 

795 

2030 

555 

461,760 

832 

APPENDIX  B 

DETERMINATION  OF  EXPECTED  TRANSIT  TIME  (ETT)  EQUATIONS 
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A.  Introduction 

In  the  Welland-Niagara  studies,  where  the  configuration  consists  of 
parallel  channels,  a  channel  choice  decision  rule  becomes  necessary  to  dis¬ 
tribute  traffic  between  the  two  branches.  The  rule  will  make  an  uncondi¬ 
tional  assignment  of  those  ships  whose  size  prohibits  them  from  using  one 
of  the  channels,  otherwise,  channel  assignment  will  be  made  on  the  basis  of 
least  expected  transit  time.  When  a  vessel  arrives  at  the  channel  choice 
point,  the  simulation  model  invokes  the  expected  transit  time  (ETT)  equa¬ 
tions  to  execute  this  rule.  The  derivation  of  the  ETT  equations  for  the 
Niagara  Canal  is  accomplished  through  the  "Experience  Data  Base"  (EDB) 
technique. 

Essentially,  the  EDB  procedure  involves  separating  a  multichannel 
system  into  distinct  single  channels.  Each  channel  is  then  simulated  in¬ 
dividually  to  determine  its  unique  operating  characteristics  under  various 
assumptions  about  traffic  levels  and  traffic  composition,  and  about 
alternative  service  facilities  and  operating  rules  within  the  branches. 

Once  a  set  of  observations  of  channel  performance  have  been  recorded, 
regression  analysis  is  used  to  relate  conditions  at  the  time  vessels 
entered  the  channel  to  their  subsequent  transit  times. 

The  EDB  procedure  is  not  used  to  develop  the  ETT  equation  for  the 
Welland  Canal,  however,  because  the  elaborate  control  system  that  current¬ 
ly  exi sts  at  Welland  augers  numerous  difficulties  for  simulation. 

Empirical  data  is  available  in  the  form  of  daily  vessel  records  and  these 
were  obtained  from  the  St.  Lawrence  Seaway  Authority. 


B.  WELLAND  CANAL  ETT  EQUATIONS 

Welland  Canal  transit  data  for  the  month  of  August,  1971  was  used  as 
the  basis  for  this  analysis.  The  format  of  this  data  as  originally  supplied 
for  downstream  travel  is  shown  in  Table  B.l.  Data  for  upstream  travel 
was  provided  in  similar  format.  Preliminary  data  treatment  consisted  of 
transforming  this  data  into  a  transit  time  matrix  as  shown  in  Table  B.2. 

This  matrix  was  then  subjected  to  regression  analysis  with  vessel  transit 
time  as  the  dependent  variable. 

For  purposes  of  regression,  it  was  found  necessary  to  model  Welland 
in  the  most  detailed  manner  possible.  To  illustrate,  on  the  macro  side, 
Welland  could  simply  be  modeled  as  one  entity,  and  the  Independent  vari¬ 
ables  in  Table  B.2  could  only  consist  of  two  rows,  viz. ,  the  number  of 
downbound  vessels  in  the  canal  and  the  number  of  upbound  vessels  in  the 
canal.  Such  a  model  was  attempted  but  the  degree  of  predictability  in  the 
resultant  regression  equation  was  extremely  small.  On  the  micro  side, 
Welland  could  be  modeled  as  the  set  of  locks,  reaches  and  bridges  as  it 
currently  exists,  and  the  independent  variables  could  consist  of  the  number 
of  upbound  and  downbound  vessels  in  each  reach  and  at  each  upbound  and 
downbound  queues  of  locks  and  bridges.  The  data  requisite  for  this  model 
was  not  available,  however.  The  most  detailed  model  that  could  be  estab¬ 
lished  consisted  of  a  six-entity  (reach)  representation  of  Welland  with  the 
first  and  last  reaches  being  the  two  waiting  areas  at  either  ends  of  the 
canal.  This  model  is  schematically  portrayed  in  Figure  2,  p.  18* 

The  conventional  linear  model  was  used  in  the  form: 


TABLE  B.l.  REPRESENTATIVE  WELLAND  TRANSIT  DATA  DAILY  RECORDS 
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where  TRj  and  represent  the  observed  transit  time  and  system  conditions 
for  vessel  j .  The  symbol  Uj  represents  a  random  disturbance  term.  The 
basic  variables  are  shown  in  the  first  column  of  Table  B.2.  These  were 
not  the  only  variables,  however,  since  power  and  logarithmic  transforma¬ 
tions  were  also  made. 

The  final  results  of  the  regression  procedure  are  shown  in  Table  B.3. 

Table  B.4  provides  a  description  of  the  variables  present  in  the  equations. 

2 

The  fraction  of  explained  variance,  R  ,  was  much  higher  when  the  dependent 
variable,  total  transit  time  (TR),  did  not  include  the  waiting  time  (see 
Table  B.l,  second  to  last  column).  Inclusion  of  waiting  time  in  the 
dependent  variable  Increased  the  degree  of  unpredictability.  This  may  be 
due  to  certain  features  of  the  existing  control  system  at  Welland,  features 
such  as  vessels  called  out  of  the  waiting  area  for  tandem  lockages,  vessel 
delays  due  to  adverse  weather  conditions,  change  of  pilot,  inspections  and 
other  aspects  of  regulated  traffic. 

C.  NIAGARA  CANAL  ETT  EQUATIONS 

EDB  simulation  runs  were  executed  for  four,  five  and  six  lock  con¬ 
figurations  with  fleet  data  generated  from  1980,  2000  and  2030  normal 
growth  projections.  The  output  from  each  run  was  again  transformed  into  a 
transit  time  matrix  such  as  that  shown  in  Table  B.2.  The  matrices  for  each 
configuration  were  merged  so  as  to  cover  the  widest  possible  range  of  sys¬ 
tem  conditions  and  were  subjected  to  regression  analysis.  The  regression 
models  for  these  Niagara  configurations  were  more  detailed  than  that  for 
Welland,  in  that  the  independent  variables  representing  system  conditions 
included  not  only  the  reach  variables  (number  of  vessels  by  direction  of 


TABLE  B. 3  WELLAND  ETT  EQUATION  (in  hours) 
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LAKE 

ERIE 


Reach  1  extends  from  Lock  1  to  Lock  7 

Reach  2  extends  from  Lock  7  to  Bridge  10 

Reach  3  extends  from  Bridge  10  to  Bridge  18 

Reach  4  extends  from  Bridge  18  to  Lock  8 


LAKE 

ONTARIO 


Length : 

vessel 

Lagnew: 

transit 

ND. : 

4 

number 

ND2: 

number 

RU6: 

square 

*V 

square 

RUX: 

square 

RD^ 

square 

SD3: 

square 

sd4j 

square 

sur 

square 

LD, : 

0 

lo*10  0 

LUX: 

lo«10  c 

lu6! 

lo®10  c 

square  of  number  of  ships  upbound  in  reach  1 
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travel)  but  also  the  lock  variables  consisting  of  the  upstream  and 
downstream  queue  sizes  at  each  individual  lock. 

The  Niagara  Canal  ETT  equations  are  shown  in  Tables  B.5,  B.6 
and  B.7  for  the  four,  five  and  six  lock  configurations  respectively. 

An  interpretation  of  the  variables  is  given  in  Table  B.8  along  with 
schematic  diagrams  of  the  relative  locations  of  locks  and  reaches 
for  each  configuration.  In  general,  the  equations  conform  to  prior 
expectations  that  the  transit  time  of  a  vessel  would  be  most  sensitive 
to  the  density  of  traffic  at  the  nearest  lock. 


TABLE  B. 5  FOUR-LOCK  NIAGARA  ETT  EQUATION  (in  minutes) 
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TABLE  B.8  NIAGARA  ETT  VARIABLES 


WELLAND 


LAKE 

ONTARIO 


WELLAND 


LAKE 

ONTARIO 


WELLAND 


LAKE 

ONTARIO 


(ill)  Six-lock  configuration  (with  nine  reaches) 


LNQi 

lfqa 

RSD 

“i 


*  number  of  vessels  in  near^  queue  at  lock  i 

-  number  of  vessels  in  far  queue  at  Lock  i 

2 

■  number  of  vessels  travelling  same  direction  in  Reach  j 

-  number  of  vessels  travelling  opposite  direction  in  Reach  j 


^ln  relation  to  the  subject  vessel 
2 

in  relation  to  the  subject  vessel's  direction  of  travel 


APPENDIX  C 


DESCRIPTION  OF  THE  QUEUING  MODEL 
USED  IN  THE  SINGLE  WELLAND  STUDIES 


i 
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The  queuing  model  used  to  supplement  the  Single  Welland  Simulation 

Studies  (see  Section  IV)  is  a  simple  waiting-line  model,  mathematically 

categorized  as  M/G/l.  That  is,  arrivals  into  the  service  facility  are 

assumed  to  be  Poisson  distributed,  there  is  only  one  service  facility  and 

the  service  times  for  the  arrival  units  are  Independent  with  some  common 

2 

probability  distribution  whose  mean  1/y  and  variance  a  are  known. 

The  terminology  used  below  is  as  follows: 

X  =  mean  arrival  rate  (expected  number  of  arrivals  per  minute) 
y  ■  mean  service  rate  (expected  number  of  units  completing  service 
per  minute) 

2 

0  »  variance  of  the  service  distribution 

p  =  A/y  ■  traffic  intensity 
W  **  waiting  time  in  queue  in  minutes 

Under  the  assumptions  given  above,  the  Pollaczek-Khintchine  formula  estab¬ 
lishes  the  waiting  time1  as: 


W 


.22  2 
A  a  +  p 

2A  (1-p) 


In  the  context  of  this  analytical  model,  the  Welland  Canal  is  viewed 
as  one  entity — a  single  service  facility.  Arrival  units  incur  delays  in 
waiting  areas  at  either  ends  outside  the  Welland  Canal  and  the  transit 
time  through  the  canal  itself  is  fairly  constant.  This  approach  was 
formulated  on  the  basis  of  the  empirical  data  as  given  below.  The  transit 
time  distribution  through  the  Welland  Canal  has  a  very  small  variance  while 
the  delay  distributions  have  large  variances. 

Although  a  state-dependent  function  for  Welland  Canal  transit  time 
was  obtained  through  regression  (i.e.,  the  transit  time  through  the  canal 


^Note:  when  the  service  distribution  is  exponential, 
the  above  equation  reduces  to  the  M/M/1  case  with 


2 

so  that  0 

X 

y(y  -  X) 


then 
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Transit  Time 

Upstream 

Welland 

Downstream 

Distribution 

Waiting  Area 

Canal 

Waiting  Area 

Mean  (min . ) 

151.00 

615.00 

149.00 

Standard 

Deviation  (min.) 

122.00 

131.00 

102.00 

St.  Dev.  as  % 

of  the  Mean 

80.79% 

21.30% 

68.45% 

Source:  VESSEL  TRANSIT 

ANALYSIS,  St.  : 

Lawrence  Seaway  Authority 

• 

was  expressed  as  a  function  of  the  system  conditions),  the  major  contri¬ 
bution  to  the  magnitude  of  transit  time  was  the  intercept  term  with  other 
variables  providing  for  minor  fluctuations  Indicating  again  that  once  a  ship 
enters  the  canal,  its  actual  transit  time  is  almost  constant. 

The  primary  objective  of  this  effort  then  became  the  prediction  of 
delays  in  waiting  areas  outside  the  canal.  The  rate  at  which  the  canal 
could  absorb  the  arrival  flows  and  dissipate  queues  needed  to  be  identified. 
This  rate,  or  more  accurately  this  service  distribution  was  determined 
through  analysis  of  the  VESSEL  TRANSIT  ANALYSIS  daily  reports  for  the  months 
of  June  and  August. 

The  traffic  flows  for  these  months  are  shown  in  Figure  C.l.  Weekly 

patterns  may  be  readily  seen,  although  there  are  enough  variations  to 

make  accurate  prediction  difficult.  Average  flows  were  25.1  and  24.7  per  day 

2 

discounting  multiple  lockages.  Since  the  number  of  daily  arrivals  vary 
greatly  from  day  to  day,  the  delay  times  for  these  units  were  associated 
with  the  traffic  density  into  the  canal  during  the  24  hours  preceding 
an  arrival. 


2 

Approximately  11%  of  the  arrival  units  were  assumed  to  undergo  multiple 
lockages . 


Figure  C.l.  Traffic  Flows  at  the  Welland  Canal 


The  results  from  this  analysis  are  shown  in  Figure  C.2.  The  service 
distribution  finally  obtained  is  described  in  terms  of  an  average  service 

O  _ 1 

cycle  ,  V  equal  to  .0112  min  (1/U  ^90  min.)  and  a  standard  deviation, 
a  equal  to  .0288  min”1.  The  waiting  time  characterized  by  these  parameters 
corresponds  well  with  empirical  data  over  a  wide  range  of  traffic. 


^ A  service  cycle  is  the  time  to  process  two  vessels  from  opposite  directions 
consecutively. 


AVERAGE  WAITi 


82 


Figure  C.2.  Comparison  of  Analytical  Results  with  Empirical  Data 
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TABLE  D.l.  WELLAND-FOUR  LOCK  NIAGARA  SUBSYSTEM 


1.  Reach  Data 

for  Niagara 

Canal 

Reach  Identification 

Length 

Mean 

Transit  Time 

Standard 

Number 

(miles) 

(minutes ) 

Deviation 

2500 

8.1 

29 

2550 

1.5 

10 

Based  on  an  analysis 

2600 

L  1 

82 

of  the  reach  between 

the  Eisenhower  &  Snell 

2650 

1.2 

24 

locks,  the  standard 

2700 

18.8 

121 

deviation  was  set  equal 
to  25%  of  the  mean.  A 

2750 

4.0 

40 

normal  distribution  was 

2800 

62.7 

221 

used. 

2.  Reaches  2100  and  2450  are  dummy  reaches  for  simulation  purposes  only. 

3.  Transit  times  through  the  Welland  reaches  were  determined  from  the 
system  conditions  during  simulation  by  way  of  the  ETT  equations. 
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TABLE  D.2.  WELLAND- FIVE  LOCK  NIAGARA  SUBSYSTEM 


LAKE 

ONTARIO 


1.  Reach  Data  For  Niagara  Canal 


Reach  ID 

Length 

Mean  Transit  Time 

Standard 

Number 

(miles  ) 

(minutes ) 

Deviation 

2500 

8.1 

29 

2550 

1.5 

10 

Set  to  252  of  the  mean 

2600 

1.7 

34 

with  a  normal  distribution. 

2650 

2.1 

43 

2700 

.4 

8 

2750 

18.8 

121 

2800 

4.0 

40 

2850 

62.7 

221 

2 .  Reaches 

2100  and 

2450  are  dummy  reaches 

for  simulation  purposes  only. 

3.  Transit  times  through  the  Welland  reaches  were  determined  from  the 
system  conditions  during  simulation  by  way  of  the  ETT  equation. 


* 
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TABLE  D.3.  WELLAND-SIX  LOCK  NIAGARA  SUBSYSTEM 


LAKE 

ONTARIO 


1 .  Reach 

Data  For  Niagara  Canal 

Reach  ID 

Length 

Mean  Transit  Time 

Standard 

Number 

(miles) 

(minutes) 

Deviation 

2500 

8.1 

29 

2550 

1.2 

8 

2600 

1.8 

36 

Set  to  252  of  the  mean  with 

2650 

0.4 

8 

a  normal  distribution. 

2700 

0.7 

10 

2750 

0.4 

8 

2800 

18.7 

121 

2850 

4.0 

40 

2900 

62.7 

221 

2 .  Reaches 

2100  and 

2450  are  dummy  reaches 

for  simulation  purposes 

only. 

3.  Transit  times  through  the  Welland  reaches  were  determined  from  the 
system  conditions  during  simulation  by  way  of  the  ETT  equations. 
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TABLE  D.4.  WELLAND  SUPERLOCKS  SUBSYSTEM 


2100 


1100 


<<- 


2200 


<  <—  <  < 


'2400  /  /2500  ,  /  2600 
-<  < - <  < - 


1200 


1300 


1400 


1500 


LAKE 

ONTARIO 


Reach  ID 

Length 

Mean  Transit  Time 

Standard 

Number 

(miles) 

(minutes ) 

Deviation 

2100 

66.7 

261 

2200 

18.8 

121 

Set  to  25%  of  the  mean  with 

2300 

.4 

8 

a  normal  distribution. 

2400 

2.1 

43 

2500 

1.7 

34 

2600 

9.6 

39 

88 


Tables  D.5  through  D.12  summarize  fleet  and  lock  data  observed  at  Poe 
Lock  (Sault  Ste  Marie)  during  the  period  May  -  July  1972. 
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TABLE  D.ll.  VESSEL  TRANSIT  BY  LENGTH  0*  VESSEL 


CL.  1  VES 

3 

0,007 

CL,  2  VES 

235 

0.503 

CL,  3  VES 

13* 

0.337 

CL. *  VES 

11 

0.027 

CL, J  VES 

1* 

0.S45 

TOTAL 

403 

1.000 

TABLE  D.12.  AVERAGE  TONS  PER  LOADED  VESSEL  BY  COMMODITY 


TONS 

0E4  LOADED 

1  VESSEL 

VESSEL 

CLASS 

OOAlT 

COAL 

PETSOLc'11  CCm'NT 

ISON  0»E 

CHE9 

HISC 

CL.l  VES 

0. 

0, 

001», 

5*00. 

0. 

0. 

4300. 

CL. 2  VES 

1463* • 

12S9S, 

7907  • 

13159, 

2095 , 

0. 

10000. 

Cl. 3  VES 

26040, 

21000. 

.  0. 

27239 , 

0. 

0. 

CL. 4  VES 

0. 

0. 

0. 

0. 

39105, 

0. 

0, 

Cl. 3  VES 

0, 

0. 

0. 

0, 

57617, 

0. 

0. 

■? 


Tables  D.13  through  D.17  contain  £leet  data 
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October  Total  904  974  1,055  1,157  1,269  1,412  1,574 

Up  462  498  539  591  648  721  804 

Down  442  476  516  566  621  691  770 
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